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ABSTRACT 

Absorption or emission against the cosmic microwave background radiation (CMB) may be observed in 
the redshifted 21cm line if the spin temperature of the neutral intergalactic medium prior to reionization 
differs from the CMB temperature. This so-called 21cm tomography should reveal important information 
on the physical state of the intergalactic medium at high redshifts. The fluctuations in the redshifted 
21 cm, due to gas density inhomogeneities at early times, should be observed at meter wavelengths by 
the next generation radio telescopes such as the proposed Square Kilometer Array (SKA). Here we show 
that the extra-galactic radio sources provide a serious contamination to the brightness temperature 
fluctuations expected in the redshifted 21 cm emission from the IGM at high redshifts. Unless the 
radio source population cuts off at flux levels above the planned sensitivity of SKA, its clustering noise 
component will dominate the angular fluctuations in the 21 cm signal. The integrated foreground signal 
is smooth in frequency space and it should nonetheless be possible to identify the sharp spectral feature 
arising from the non-uniformities in the neutral hydrogen density during the epoch when the first UV 
sources reionize the intergalactic medium. 

Subject headings: (cosmology:) early universe — galaxies: general 



1. INTRODUCTION 

Studies of Lya forest in the spectra of high redshift 
QSOs have now clearly established that the Universe is 
highly ionized by z ~ 6 and that the intergalactic medium 
(IGM) developed extensive nonlinear structures at these 
high redshifts. Although it is now possible to simulate 
the formation of the first ionizing objects in the universe 
(Gnedin & Ostriker 1997, Abel et al. 1998, Abel, Bryan 
& Norman 2000, Gnedin 2000, Gnedin & Abel 2001), the 
epoch at which these objects were formed, their nature and 
the timescale over which the transition from the neutral 
to the reionized universe occurs is currently unknown. 

Building on the pioneering work of Field (1958, 1959) 
Scott and Rees (1990), Madau, Meiksin and Rees (1997; 
MMR97 hereafter) proposed that 21 cm tomography can 
probe the IGM prior to the epoch of full reionization 
(z > 6). The idea is that the radiation from the first dis- 
crete sources should be detected indirectly through its in- 
teraction with the surrounding neutral IGM and in the re- 
sulting emission and/or absorption against the cosmic mi- 
crowave background (CMB) at the frequency correspond- 
ing to the redshifted 21 cm line. The resulting patchwork 
(both in angle and in frequency) in the 21 cm radiation 
resulting from non-uniformities in the distribution of gas 
density and distribution of Lya sources (i.e., the "Cosmic 
web" at these early times) should be measurable with the 
next generation of radio telescopes (e.g., Tozzi et al. 2000; 
hereafter TOO) such as the proposed Square Kilometer Ar- 
ray (SKA) 5 or the Low Frequency Array (LOFAR) 6 . 

In this paper, we assess the extent to which the detection 



of the redshifted 21 cm emission fluctuations is impeded 
by confusion noise from extra-galactic foreground of radio 
sources. 

In the following section we briefly summarize the phys- 
ical origin and expected magnitude of the 21cm signa- 
ture. Section ^ describes the extra-galactic radio fore- 
ground emission and calculate its expected confusion noise 
component. Finally, in §4 and §5 we present and discuss 
our results. 

2. THE EXPECTED SIGNAL 

Detailed investigations of the absorption/emission sig- 
nal expected in the 21 cm radiation have been carried out 
by MMR97 and TOO. It is anticipated that a strong soft 
UV background radiation field is established by the first 
stars and quasars prior to reionization (e.g., Haiman, Abel, 
and Rees 2000). Lyman alpha photons from this back- 
ground will couple the spin temperature of the neutral 
IGM to its kinetic gas temperature predominantly via the 
Wouthuysen-Field effect (Wouthuysen 1952; Field 1958). 
Since the IGM thermally decouples from the CMB at red- 
shift ~ 130 and cools adiabatically thereafter (Tigm ^ 
(1 + z) 2 ) one finds Tjgm/Tcmb "C 1 before the first lu- 
minous objects reheat the universe. Here the signal is in 
absorption against the CMB, and it is expected at the 
wavelength of 21cto x (1 + ZL ya ), where zj, ya denotes the 
redshift at which the spin temperature is first coupled to 
the kinetic gas temperature. Once the IGM is heated to 
a temperature above that of the CMB, the signal will be 
in emission and roughly independent of the spin tempera- 
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ture. The differential antenna temperature (the brightness 
temperature T B = T CMB e~ T + T S [1 — e~ T ] ) observed at 
Earth between such a patch of IGM at a spin temperature 
Ts and the CMB is approximately given by (TOO), 

where h denotes the current Hubble constant in units of 
100 km s -1 Mpc -1 , and Ob gives the baryon density in 
units of the critical density. \ST\ is larger by a factor 
Tcmb/Ts (^18 for reionization redshifts > 6) for the ab- 
sorption signal than for emission. Unfortunately, a rel- 
atively large radiation flux in Lya is required to ensure 
the coupling of the spin temperature to the kinetic gas 
temperature. Therefore, the recoil from Lya scatterings 
most likely heats the IGM to or above the CMB tempera- 
ture before a sufficient Lya background flux is established 
(MMR97). Consequently the largest signals predicted by 
TOO and MMR97 are <10 mK (at e.g., 150 MHz if z ~ 9), 
as expected from equation (|l|) and the "Cosmic Web" at 
these times is most likely to be probed in emission at 21 
cm. The appropriate resolution for measuring fluctua- 
tions in the IGM is of the order of 1 arc-minute with the 
frequency window of 1 MHz around 150 MHz (MMR97; 
TOO). However, if the absorption feature were observed, 
it would give valuable insights to the epoch at which the 
first stars formed (MMR97; TOO). At the epoch of reion- 
ization "breakthrough" the 21cm signal from diffuse gas 
decreases steeply. One should be able to observe this sig- 
nal integrated over most of the sky even with moderate 
instruments (Shaver et al. 1999). 

In the following sections we examine the confusion noise 
introduced at 150 MHz by extra-galactic radio sources and 
discuss how this might impose serious limitations for the 
21-cm tomography. 

3. EXTRA-GALACTIC RADIO FOREGROUNDS 

SKA is planned to operate in the frequency range 0.01- 
20 GHz with an angular resolution of about 1-10 arc- 
seconds and a sensitivity down to a few tens of nJy Lim- 
itations on the sensitivity achievable for the measurement 
of redshifted 21-cm signal will be set by the contamina- 
tion from the galactic and extra-galactic foregrounds. The 
dominant galactic contribution is the synchrotron back- 
ground which comprises a fraction of about 70% at 150 
MHz. On the angular scales that are most relevant for 
the 21cm tomography, we nonetheless expect the domi- 
nant angular fluctuations to be caused by clustering and 
discreteness of extra-galactic sources (discussed below) . 

The galactic foreground could only make a compara- 
ble contribution on sub-arc-minute scales if there were, 
throughout the interstellar medium and galactic halo, fluc- 
tuations of order unity in its volume emissivity. Such fluc- 
tuations would exist, but only in small regions (e.g. ac- 
tive supernova remnants). On larger angular scales, where 
the galactic foreground fluctuations would be relatively 
more important, it has been sufficiently well-studied (al- 
though at much higher radio frequencies than considered 
here) because of its importance for the observations of the 



cosmic microwave background radiation (see references in 
Tegmark & Efstathiou 1996). A detailed discussion of the 
Galactic spectral contamination to the redshifted 21-cm is 
given in Shaver et al. 1999. 

Detailed multi-wavelength observations of the galactic 
radio emission could be modeled sufficiently accurately for 
our purpose (at least in some regions of the sky where ob- 
servations would then take place). Hence, we limit our dis- 
cussion to the confusion noise introduced by extra-galactic 
foreground sources such as radio galaxies, AGN and nor- 
mal galaxies which are likely to dominates the radio counts 
at the low flux density levels. 

3.1. Counts of radio sources at low radio frequencies 

The wavelength region of interest is from 50 MHz (z ~ 
30) to 200 MHz (z ~ 6). To evaluate the impact of extra- 
galactic foreground radio sources in this wavelength range 
it is necessary to model the number density of sources as 
a function of flux (the differential counts N(S) per stera- 
dian). At present, the appearance of the radio sky at flux 
density levels below 1/z Jy is not well known. However 
deep VLA surveys have allowed to extend direct deter- 
minations of radio source counts down to a few /xJy (at 
v > 1.4 GHz), implying a coverage of about 7 orders of 
magnitude in flux. At lower frequencies the limiting flux 
densities are higher because of the confusion noise intro- 
duced by extended sources. 

The source counts from the 6C survey (Hales, Baldwin 
& Warner 1988) which was carried out at 151 MHz are a 
useful guide. The radio source counts at all frequencies are 
typically well described by a Euclidean power law region 
at the highest flux densities followed by a flatter portion at 
lower flux densities (e.g. Formalont et al. 1991). We there- 
fore extrapolate the 151 MHz differential source counts by 
a double power-law fitted to the observed counts. This 
gives: 

where 7l = 1.75, 72 = 2.51, h = k 2 S^ 2 ~ J1 with k 2 = 4.0 
per sr per mJy and So = 880 mJy. 

This fit also includes the counts from the 3CR survey 
and the 3 CRR catalogue at 178 MHz (Laing, Riley & 
Longair 1983) transposed to 151 MHz assuming a mean 
spectral index a = 0.75 (S oc v~ a ), typical for the emis- 
sion from extended lobes at these frequencies (e.g.; Laing, 
Riley & Longair 1983). The limiting flux density for these 
survey was ~ 100 mJy and our extrapolation is somewhat 
uncertain. However, we will show that our main conclu- 
sions are insensitive to the particular choice of the source 
counts. 

3.2. Sky fluctuations from unresolved sources 

The contribution to CMB fluctuations from randomly 
distributed sources of various nature has been extensively 
discussed in the literature (e.g. Scheuer 1974 Cavaliere 
& Setti 1976; Franceschini, Vercellone & Fabian 1998; 
Tegmark & Efstathiou 1996; Scott & White 1998; Tof- 
folatti et al. 1999; Perna & Di Matteo 2000). We can ap- 
ply these techniques directly to the problem at hand. We 
estimate the confusion noise due to unidentified sources 
below a flux density limit S cu t which will contribute to 
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fluctuations in the energy band of redshifted 21-cm emis- 
sion. In the absence of clustering (Poisson distribution), 
the angular power spectrum, C/, is a white noise (C; ~ 
constant). The contribution to the background below the 
flux cut S cu t due to sources with a Poisson distribution is 
given by: 
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However, the analysis of large samples of nearby radio- 
galaxies has shown that sources are typically strongly clus- 
tered (e.g. Peacock & Nicholson 1991). Clustering de- 
creases the effective number of objects in randomly dis- 
tributed pixels and consequently enhances the pixel-to 
pixel fluctuations (e.g.; Peebles 1980). In the case of a 
powerlaw angular correlation function (w(8) — (9/9q)~P ) 
the power spectrum of intensity fluctuations due to clus- 
tered sources can be simply estimated as (e.g. Scott & 
White 1999) 

C cluster = W[I 2 ^ 

where wi cx l^~ 2 is the Legendre transform of w(9) and 
I = L c "* S(dN/dS)dS is the background contributed by 
sources below S cu t- If sources are clustered like galaxies 
today or as Lyman-break galaxies (Giavalisco et al. 1998) 
at z ~ 3 we expect (3 w 0.8 — 0.9 (most of the sources at 
the faint flux density levels of interest here would in-fact 
be galaxies). The angular correlation 9q, is roughly in- 
dependent of redshift because of the bias of high redshift 
halos and typically of the order of a few arc-minutes (see 
e.g., Oh 1999). Here we take 9$ = 4 arc-minutes. 

We can compute the angular rms temperature fluctua- 
tion from 
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is the conversion factor from temperature to flux density 
(per steradian), B V (T) is the Planck function, x = hv/kT 
with T = 2.725 K (Mather et al. 1999) the CMB tempera- 
ture. Note that the temperature fluctuations increase with 
increasing angular resolution. 

4. CONFUSION NOISE 

The rms brightness temperature fluctuation for various 
values of S cu t, as a function of angular scale 9, is shown in 
Figure 1. The thicker lines represent the clustering term 
and the thinner lines the Poisson term. 

The foreground signal from the point sources is typi- 
cally larger than the lOmK signal (eq. [TJ) expected from 
the 21 cm radiation from the IGM at high redshifts. As 
expected, the signal decreases for decreasing values of S cut , 
i.e. if more sources can be identified and removed. In par- 
ticular, as the value of the cutoff flux S cu t is lowered, the 
Poisson term declines a lot more sharply than the clus- 
tering term, until the latter totally dominates. Given our 
dN/dS, both components are dominated by objects just 



below the detection threshold S cu t , but this dependence is 
stronger for the Poisson term. Thus subtraction of bright 
sources decreases the Poisson term much more effectively 
than the clustering term. 

The SKA rms sensitivity is estimated to be (Taylor & 
Brown 1999) 7 
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at 150 MHz. Therefore, in flux density units, the in- 
strumental noise Sinstr is independent of angular resolu- 
tion. Given the appropriate angular resolution required 
to beat down source confusion one can therefore identify 
and remove foreground sources down to flux densities of 
e.g. S'cut.iim = 7Sinstr- A limit on the angular resolution 
required to identify all sources down to S C ut,Hm can be de- 
termined by comparing S C ut,Hm to the confusion noise due 
to sources below this cutoff: 

c c / o \ °- 75 / a \ -1-55 

«J <->cut -, n / '-'cut \ | 

N ~ {Cf luster \\ s t 2 )V2 ~ l,2.1//Jy J \073" 

(8) 

Given that SKA is planned to achieve an angular resolu- 
tion of the order of a fraction of arc-second, in principle it 
may be possible to remove all sources down to S ~ Scut.iim 
without being limited by the confusion noise of the sources 
below that flux density. Figure 1 shows that even if sources 
are subtracted down S C ut = 1/x Jy ~ Scutum, the fore- 
ground component (the Poisson one for 9 < 10 arc-minutes 
and clustering one for any 9) still dominates the primary 
21 signal. 

Note also that S'cut.iim will not be achieved if sources are 
even slightly extended, i.e. if their angular extent 9 > 0.3 
arc-seconds which is a small number even at large dis- 
tances. At these faint flux density levels intrinsic source 
confusion due to their finite sizes would become important 
(see e.g., Kellermann & Richards 1999) 

These results clearly depend on the extrapolation of the 
radio source population (Eq. 2) down to flux density levels 
of the order of SKA sensitivity. Clearly, if the radio source 
population were to cut-off at flux densities above the SKA 
detection threshold, then all sources could be identified 
and removed leaving no foreground contamination to the 
21-cm signal. 

However, this is unlikely given the typical shallow spec- 
tral slopes of radio sources and in particular the results 
from the deep VLA surveys, which found that radio counts 
at higher frequencies extend all the way down to the cur- 
rent instrumental sensitivities (e.g. at 5 GHz to flux den- 
sities ~ 1/^Jy) without any evidence of a turn over. In- 
deed, VLA observations have revealed a further upturn of 
the differential counts of compact sources below a few /xJy 
(Mitchell & Condon 1985; Windhorst et al. 1985) imply- 
ing A^(> S) cx S- 1 - 2 . 

5. USING FREQUENCY INFORMATION 

As we have seen in the previous section, within a re- 
alistic beam size one is integrating over a very large 
number of unresolved extra-galactic sources. Their in- 
dividual spectra are rather shallow and typical ob- 
served flux densities of galaxies, supernovae remnants 
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and the sky background etc. scale with v 13 with 
0.8 < 13 < 0.2 (Zombcck 1982). Their inte- 
grated flux will hence also have a rather shallow slope. 
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Fig. 1. — Angular dependence of the temperature fluctuations 
due to foreground sources at 150 MHz. Here the signal is shown 
after removal of all sources with flux densities S > S C ut as indi- 
cated in the figure. For each value of S cu t, the thin lines show the 
Poisson noise, while the thick lines show the corresponding noise if 
the sources are clustered. For Scut $0.1 mJy, the clustering term 
dominates the Poisson term at the scales of interest for the 21cm 
observations. Note that the largest expected signal of 21cm emission 
at high rcdshifts is below lOmK indicated by hatched region. 

The 21 cm contribution, however, has steep spectral 
structures. The spectral index along a line of sight changes 
by an amount of order unity over a small frequency range 
Ai/, where t\vjv is the scale of neutral hydrogen non- 
uniformities relative to the Hubble radius. This change 
may occur multiple times along some lines of sight because 
of ionized regions prior to ionization "breakthrough" . Sub- 
tracting out the measured powerlaw of the foregrounds 
one should see the left-over frequency dependence of the 
21 cm contribution. Since the 21 cm signal is so sharp in 
frequency space, small errors in the measured integrated 
foreground spectrum will not hinder this measurement. 
For example, observing with 2 MHz resolution with an 
error of A/3 ~ .05 of the measured foreground spectrum 
should allow to distinguish a signal that is (v/uq) ai:> — 1 = 
(152/150) A/3 - 1 - 6 x 10~ 4 times smaller than the fore- 
ground. Note that this argument assumes that the beam 
shape does not change for observations at a different fre- 
quency. Obviously, if the beam changes just slightly in size 
the clustering of the unresolved foregrounds will introduce 
noise in the spectra on the order of the change of area 
of the beam. However, it is not clear whether the large 
arrays needed to detect the signal will allow such a good 
control over the synthesized beam and a more detailed 
study of these technical limitations seems warranted. If 



these issues could be overcome, this signal would give very 
strong constraints on the nature of the sources rcionizing 
the intergalactic medium. Not only could one pin down 
the epoch of reionization but one could map the typical 
size of the regions of influence of individual sources. 

6. CONCLUSIONS 

Radio observations at meter wavelengths with instru- 
ments such as the SKA should map the large scale struc- 
tures at high redshift when the IGM had no yet been 
reionized (6<z<30). They should therefore provide a use- 
ful tool for probing the epoch, nature and sources of the 
reionization in the universe and their implications for cos- 
mology. 

Assuming that the radio Galactic foreground can be suf- 
ficiently modeled out (see Shaver et al 1999, for possible 
strategies to do so), we have shown that the confusion 
noise provided by extra-galactic radio sources provides a 
serious contamination to the brightness fluctuations ex- 
pected in the redshifted 21 cm emission from the IGM at 
high rcdshifts. In particular, even if the radio source pop- 
ulation is fully identified and removed down to the flux 
densities corresponding to the planned sensitivity of SKA, 
its clustering noise component will totally dominate the 21 
cm signal at all scales and its Poisson component at scales 
9 < 10 arc-seconds. Therefore, the detection of brightness 
fluctuations in the redshifted 21-cm seems unfeasible. 

Note however that even in the case of significant fore- 
ground brightness fluctuations the rise and decay of the 
21 cm emission is expected to show very sharp features in 
frequency space which are potentially measurable even in 
the presence of the strong extra-galactic (or galactic) con- 
tinuum. The integrated spectrum of extra-galactic fore- 
ground sources will have a featureless powerlaw energy 
spectrum which will be measured directly and can mod- 
eled out in frequency space. To detect the spectral devia- 
tions from the redshifted 21 cm signal one would require 
observations in a sufficiently narrow bandwidth <5MHz. 
In particular, given such a narrow bandwidth, searching 
for the 21-cm signal should be possible even if there are 
spectral variations in the foregrounds. Any variation in 
frequency space in the foregrounds would be correlated 
(as the galactic or extragalactic foreground have power 
law spectra) whereas the 21 cm feature would appear as 
an uncorrelated signal and should therefore be detectable. 
Technical challenges as e.g. techniques that minimize the 
change in the effective beam size as one changes the fre- 
quency band will have to be addressed. 
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